Light bullets in quadratic media with normal dispersion at the second harmonic 
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Stable two- and three-dimensional spatiotemporal solitons (STSs) in second-harmonic-generating 
media are found in the case of normal dispersion at the second harmonic (SH). This result, surprising 
from the theoretical viewpoint, opens a way for experimental realization of STSs. An analytical 
estimate for the existence of STSs is derived, and full results, including a complete stability diagram, 
are obtained in a numerical form. STSs withstand not only the normal SH dispersion, but also 
finite walk-off between the harmonics, and readily self-trap from a Gaussian pulse launched at the 
fundamental frequency. 

PACS numbers: 42.65.-k 42.65.Tg 



A multidimensional soliton is a self-supporting solitary 
pulse resulting from the balance between nonlinearity of 
the medium and diffraction and dispersion of the wave 
field. The most physically relevant realization of this is 
provided by optical spatiotemporal solitons (STSs), alias 
"light bullets" , which are self-confined in the longitudinal 
and transverse directions 0. 

Spatial optical solitons, which are not localized in the 
longitudinal direction, have been studied in depth Q, 
and observed in both one- and two- 0] dimensional 
(ID and 2D) geometries. Quasi-2D STSs, which are lo- 
calized in the longitudinal and one transverse directions, 
have been observed too [fjj. A challenge for the experi- 
mentalist is creating a fully 3D bullet. The issue has a 
special physical purport, as, unlike 2D solitons which oc- 
cur in various media, optical crystals provide for a single 
real possibility to create 3D solitons in classical physics. 
Besides their significance as fundamental objects, STSs 
may provide THz switching rates in optic-logic systems, 
if used as information bits |fj. 

It is well known that STSs are unstable in Kerr me- 
dia • They can be stabilized if the nonlinearity is sat- 
urable || . Other effects can also stabilize a spatiotempo- 
ral pulse, as shown in an experiment demonstrating sta- 
ble space-time focusing in a planar waveguide Q. How- 
ever, it was not a soliton, as the stabilization was due to 
multi-photon absorption and intra-pulse Raman scatter- 
ing, which are dissipative effects. 

A promising path for making stable STSs is to em- 
ploy second- harmonic generation (SHG). Various types 
of solitons in SHG materials have attracted a lot of in- 
terest [Io|. In particular, it was proven long ago 0, by 
means of variational estimates, that SHG must give rise 
to stable 2D and 3D solitons realizing a minimum of the 
Hamiltonian. 

Most theoretical works on STSs in SHG media assumed 
spatiotemporal isotropy 0], which implies equal GVD 
(group-velocity-dispersion) coefficients at the fundamen- 
tal frequency (FF) and second harmonic (SH), that does 



not take place in reality. A general non-symmetric case 
was considered in Ref. Using the variational ap- 

proximation (VA), STSs were predicted, and numeri- 
cal simulations verified their stable existence in the 2D 
case. Later, numerically exact 2D and 3D soliton solu- 
tions were found and their stability was confirmed |l4| . 
Then, another ingredient important for the description 
of the real physical situation was introduced, namely, the 
group-velocity mismatch (GVM, alias walkoff) between 
the harmonics, showing that stable STSs exist in this 
case too 

It is obvious that solitons in SHG media may only exist 
if GVD is anomalous at FF. Even without postulating the 
above-mentioned spatiotemporal isotropy, all the previ- 
ous works assumed that GVD must be anomalous at SH 
too. However, this assumption ignores a fundamental 
problem: in available SHG materials, FF suffers strong 
absorption if the SH wavelength is long enough for the 
dispersion to be anomalous. Note that, in Ref. VA 
formally predicted STSs for normal GVD at SH, but sim- 
ulations led to conclusion that it could really exist only 
with anomalous GVD at SH, while the prediction for the 
normal SH dispersion was an artifact of VA based on the 
Gaussian ansatz, which incorrectly treats exponentially 
decaying tails of the soliton (see a discussion of the role of 
the tails below). In the 3D case, VA predicted that STSs 
could exist at small normal values of GVD at SH, but 
this was never checked. Besides its crucial importance 
for the creation of LBs in experiment, the issue is also a 
challenge for analysis of multidimensional nonlinear-wave 
models. 

In this work we demonstrate that, contrary to the com- 
mon belief, STSs in the 3D and 2D cases effectively exists 
and can be stable in SHG media with normal GVD at 
SH. We also conclude that it survives in the presence of 
GVM. These results greatly enhance the likelihood that 
STSs will be generated experimentally. 

The co-propagation of the FF and SH waves in SHG 
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media is described by known equations 0, 0] , 

iv,£ + + u TT — u + u*v = 0, (1) 
+ W 2 j_v + 5v TT - i(iv T - jv + (1/2) u 2 = 0. (2) 

Here, u = V2£ ie - tz / z ° z lu 2 x (2) /(he 2 ) and v = 
£ 2 e~ 2lZ/ ' z ° Zotjfx I '(kic 2 ), E\,2 are complex electric-field 
envelopes at the frequencies w\ and 2oj\, ki = k (u>i) 
is the FF carrier wave number, £ = z/zo and r 
: / — z/v g )/to are the normalized propagation distance 
and reduced time, v g and to being the group velocity 
and time scale at FF, zq — 2t 2 t /k 1 , where k l is the 
GVD coefficient at FF, and the transverse radial coor- 
dinate is r — \j2k\ (x 2 
ator, Vi = d 2 /dr 2 



y 2 ) / zq. The diffraction oper- 
[(D - 2)/r)d/dr, where D = 3 
or 2 is the spatial dimension, acts on the transverse co- 
ordinates {x,y} = y/2k\) 'zqt. In the experiment, units 
of time and transverse and propagation distance in the 
scaled equations typically correspond to 100 fs, 50 fim, 
and I cm, respectively |5J. Further, a in Eq. J2J) is 
GVM, and 7 = 4 + 2z [2fc(wi) - k(2io x )\ is the phase 
mismatch between the harmonics. The equations con- 
serve the Hamiltonian, momentum, axial angular mo- 
mentum, and the Manley-Rowe invariant (energy) , which 
is Q = 27r/ °° rdr dr (\u\ 2 + 4\v\ 2 ) in the 3D case. 

In Eq. 0J, the FH dispersion is anomalous, while the 
ratio k 2 /k 1 = 5 of the SH and FH GVD coefficients is 
assumed negative, corresponding to normal GVD at SH. 
This case, which corresponds to the experimental reality, 
has not yet been studied (except for VA, with a = 0, in 
Ref. 0). 
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FIG. 1: The energy of stationary spatiotemporal solitons vs. 
the phase mismatch 7 in the case of normal GVD at the sec- 
ond harmonic (S — —0.01). The curves with 5=1 represent 
earlier known solutions and are included for comparison. 

First, an estimate for the existence of STSs can be ob- 
tained for small \S\. To this end, Eq. @ is treated as 
a linear equation for SH, driven by the term (1/2) u 2 . 
A SH field component of the CW (continuous-wave) 
form, B cw exp (ik • r — itor) , which does not vanish as 



r, \t\ — > 00 and thus prevents the existence of the soli- 
ton, is determined by zeros of the denominator of the 
corresponding Green's function, which are located at 



\S\uj 2 - auj - k 2 - 7 = 0. 



(3) 



On the other hand, an asymptotic form of the solution 
to Eq. Q at r, |r| — * 00 is u ~ p^( D_1 )/ 2 exp (— p), with 

P ■ 



= Vt 2 + r 2 . 




FIG. 2: Stability regions for 2D and 3D solitons without 
walkoff (a = 0). The solitons are stable above the respective 
borders. The insets show central profiles of the 3D soliton 
in the temporal and spatial directions, for 8 = —0.02 and 
7 = 5. The upper and lower curves pertain to the FF and SH 
components u and v, respectively. 

An estimate for the maximum amplitude B cw of the 
possible CW component can be derived. To simplify the 
presentation, one can adopt a model form of the FF field 
complying with the above asymptotics and securing a 
globally correct shape of the solution. In the 3D case, 
it is u = const • p~ 1 (sinhp) sech 2 p, but the final expo- 
nential estimate (see below) does not predicate on this 
"ansatz". Straightforward calculations [Fourier trans- 
form of u 2 (r, r), inverse transform of its product with 
the Green's function of Eq. J5J, and isolating a CW con- 
tribution from uj determined by Eq. 0] show that, in 
the case a = (no- walkoff ) , the CW intensity is expo- 
nentially small (for small \S\) if 7 > 0, an upper estimate 
for it being 



B£ w < const 



7(7/l<5|) 3/2 exp(^V7/^)- (4) 



While the pre-exponential part of this estimate 
may be ansatz-dependent, the crucially important 
exponential part is not, and it applies to the 
2D case as well; if a 7^ 0, it is replaced 

i.e., the 



by exp < — 7r 



( 7 /|*|) + (a/2Sy-\a/2S\ 



walkoff makes the estimate weaker. 

Thus, although it is not guaranteed that the STS exists 
as a rigorous solitary-wave solution, the one with the 
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CW tail obeying the exponentially small estimate may be 
tantamount to a true soliton in any possible experiment. 
The actual region of the existence and stability of STSs 
should be found numerically, which is done below. 
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FIG. 3: Evolution of the 3D soliton self-trapping from an 
initial Gaussian pulse for 8 — —0.01,7 = 1, and a — 0. 
The propagation distance corresponds to ~ 800 dispersion 
lengths. The inset shows the evolution of the FF amplitude on 
a shorter scale, for the same case and the one with 5 — 0.095, 
closer to the instability border. 

Stationary solutions were found by solving the z- 
independent version of Eqs. Q and J2J) with 6 < 0. 
Although the above analysis suggests that STS may 
have a small nonvanishing tail in the SH component, no 
tails were found, with the numerical accuracy available 
[0(1(T 4 )], in all the pulses which were identified as soli- 
tons. It was also verified that the shape of the soliton 
does not change with increasing the numerical accuracy 
and/or size of the integration domain (for instance, there 
was no change with the increase of the number of points 
from 64 3 to 256 3 ). Insets in Fig. El show an example of 
a 3D soliton, found for 5 = —0.02, 7 = 5 and a = 0. 
An STS family was generated by varying S and 7, Fig. 
n showing the soliton's energy versus 7. A conclusion of 
direct relevance to experiment is that, in the most chal- 
lenging 3D case, energy required to create a soliton with 
normal GVD at SH is smaller than in the previously 
studied case 5 — 1 by a factor ~ 2.5. 

Further consideration shows that, at large 7, the bulk 
of the STS's e nerg y is in FF (in compliance with the 
cascading limit jlfj), while at smaller 7 the amplitudes of 
both components are nearly equal, hence SH has roughly 
four times the energy of the fundamental. This shows 
that, remarkably, the soliton is able to exist keeping ~ 
80% of its energy in the normal-GVD component. 

A crucially important issue is stability of the soli- 
tons. First of all, it may be judged using the Vakhitov- 
Kolokolov criterion, which states that solitons may only 
be stable if dQ/d'j > [ljj. Figure ^ shows that the soli- 
ton families satisfy this criterion, unless 7 is too small. 



A systematic stability test was performed by direct sim- 
ulations of Eqs. and J5J. The results are summarized 
in Fig. El 

At small \6\, the stability border in Fig. 121 is linear, 
which complies with the fact that the exponential factor 
in the estimate @ is a function of the ratio j/\S\ (the ex- 
ponential estimate yields a factor ~ 10~ 5 at the stability 
border in the figure). Furthermore, the 2D and 3D stabil- 
ity borders are almost identical at small |<5|, in agreement 
with the fact that the exponential estimate is the same for 
both dimensions. The solitons are truly robust: not only 
they are not destroyed by perturbations, but they also 
readily self-trap from initial pulses of quite an arbitrary 
shape, which is illustrated by Fig. |31 Note that the tran- 
sition from the initial profile to the soliton's one takes the 
propagation distance of few dispersion lengths [zd), and 
the STS remains stable over an extremely long distance 
ph 800 Zd- This shows not only that the STS will be stable 
in any experimental setup, but also that possible energy 
leakage due to the formation of the above-mentioned CW 
tail cannot be spotted in the course of the extremely long 
propagation. The same figure shows persistent intrinsic 
oscillations in the stable STS, which is attributed to ex- 
citation of an intrinsic mode. The excitation is tangible 
if the initial pulse is launched at a point which is located 
deep within the stability region. If the pulse is taken 
close to the stability border, the inset to Fig. |31 shows 
that the vibrations quickly fade, the pulse relaxing to a 
stationary soliton through transient emission of radiation 
(the same is observed in the 2D case). 
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FIG. 4: The evolution of the FF amplitude for 3D solitons 
with 7 = 15,5 = —0.1 and various values of the GVD pa- 
rameter a. The soliton is generated from a Gaussian pulse 
launched in the FF component. It exists up to a — 3. 

The realistic model must include GVM. A previous 
work addressed "walking" STSs, but in the case S > 
0. Starting with a Gaussian pulse in the FF component, 
we tried to generate STSs with a finite GVD parameter 
a. Figure 01 shows the result: at large 7 and relatively 
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large negative S, the bullet self-traps and is found to be 
fully localized for moderate values of a. The pulse cannot 
self-trap (decaying into radiation) if a > 2. The stability 
diagram displayed in Fig. |2]can be extended, adding a as 
a parameter. For a ^ 1, the stability region is nearly the 
same as at a = 0, but with the increase of a, it quickly 
shrinks and no longer exists for a > 2. 

Lastly, stationary 3D solitons with an internal vortic- 
ity were found too. However, they are always unsta- 
ble against perturbations breaking their axial symmetry, 
similar to what was found for the same type of STS in 
the case 5 > |T^ |. 

In conclusion, we have shown in the numerical form 
that 2D and 3D stable spatiotemporal solitons exist in 
quadratically nonlinear media when the dispersion is 
anomalous at FF but normal at SH. An exponential an- 
alytical estimate for the existence of the solitons was ob- 
tained. The solitons readily self-trap from a Gaussian 
pulse launched at FF, showing no trace of decay over ex- 
tremely long propagation distance. The solitons have the 
energy smaller than their counterparts in the case of the 
anomalous dispersion at SH, and they exhibit robustness 
to the walkoff between the harmonics. These features 
imply that the solitons of this type can be created under 
a variety of experimental conditions. 

We appreciate a valuable discussion with D. Mihalache 
and K. Bcckwitt. The work was supported in a part by 
the Binational (US-Israel) Science Foundation (contract 
No. 1999459) and by a matching grant from the Tel 
Aviv University (TAU). Access to supercomputing facil- 
ities was provided by the High Performance Computing 
Unit at TAU. 



* Current address: School of Mathematics and Statistics, 
University of New South Wales at ADFA 2600, ACT, 
Australia 

[1] Y. Silberberg, Opt. Lett. 15, 1281 (1990). 
[2] G.I. Stegeman and M. Segev, Science 286, 1518 (1999). 
[3] R. Schiek, Y. Baek, and G.I. Stegeman, Phys. Rev. E 53, 
1138 (1996). 

[4] W.E. Torruellas et al, Phys. Rev. Lett. 74, 5036 (1995). 
[5] X. Liu, L.J. Qian, and FW. Wise, Phys. Rev. Lett. 82, 

4631 (1999); X. Liu, K. Beckwitt, and FW. Wise, Phys. 

Rev. E 61, R4722 (2000); X. Liu, K. Beckwitt, and FW. 

Wise, Phys. Rev. Lett. 85, 1871 (2000). 
[6] R. McLeod, K. Wagner, and S. Blair, Phys. Rev. A 52, 

3254 (1995). 

[7] J.J. Rasmussen and K. Rypdal, Phys. Scripta 33, 481 
(1986). 

[8] A.B. Blagoeva et al, IEEE J. Quantum Electron. QE- 

27, 2060 (1991). 
[9] H.S. Eisenberg et al, Phys. Rev. Lett. 87, 043902 (2001). 
[10] C Etrich et al, Progr. Opt. 41, 483 (2000); A.V. Buryak 

et al, Phys. Rep. 370, 63 (2002). 
[11] A. A. Kanashov and A.M. Rubenchik, Physica D 4, 122 
(1981). 

[12] D.V. Skryabin and W.J. Firth, Opt. Commun. 148, 79 
(1998). 

[13] B.A. Malomed et al, Phys. Rev. E 56, 4725 (1997). 
[14] D. Mihalache et al, Opt. Commun. 152, 365 (1998); D. 

Mihalache et al, Opt. Commun. 159, 129 (1999). 
[15] D. Mihalache et al, Opt. Commun. 169, 341 (1999); D. 

Mihalache et al, Phys. Rev. E 62, 7340 (2000). 
[16] D. Mihalache et al, Phys. Rev. E 62, R1505 (2000). 



